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The differences of the masses of isotopes with atomic numbers between ~ 10 and ~ 30 can be 
described within the chiral soliton model in satisfactory agreement with data. The rescaling 
of the model is necessary for this purpose — decrease of the Skyrme constant by ~ 30%, 
providing the "nuclear variant" of the model. The asymmetric term in Weizsacker — Bethe — 
Bacher mass formula for nuclei can be obtained as the isospin dependent quantum correction 
\ to the nucleus energy. Some predictions of the binding energies of neutron rich isotopes are 

■ made in this way from, e.g. 16 Be, 19 B to 31 Ne or 32 Na. The neutron rich nuclides with 

high values of isospin are unstable relative to the decay due to strong interactions. The SK4 
Cln 1 (Skyrme) variant of the model, as well as SK6 variant (sixth order term in chiral derivatives 

in the Lagrangian as solitons stabilizer) are considered, the rational map approximation is 
used to describe multiskyrmions. 
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1. Introduction 

o 

In the absence of the complete theory of strong interactions and nuclear forces the 
checking of fundamental principles, which are believed to hold in any theoretical model, 
can be useful and of great importance. The description of some well established and rather 
general properties of nuclei, as well as searches and studies of unusual forms of matter, in 
particular the neutron rich nuclides, may provide important source of lacking information 
and impact for the development of new concepts and ideas. 

The effective field theories are the powerful tool for the studies not only of mesons, 
baryons and their interactions at low energy [1—3], but also baryonic systems (nuclei) which 
appear as quantized bound states of skyrmions with baryon number B > 2. The properties of 
the deuteron and 3 He, 3 H and 4 He were explained semiquantitatively [4—5] starting from few 
basic principles and ingredients contained in the effective chiral Lagrangian of the Skyrme 
model [1] or its modifications (short review of early results can be found, e.g. in [6]). 

For the B = 2 system the quantization of the bound state of skyrmions possessing 
originally characteristic torus-like form [7] allows to obtain the deuteron and the singlet NN- 
scattering state with isospin / = 1 [8] . The binding energy of the deuteron is about ~ 30 MeV 
when zero modes quantum corrections are taken into account. It decreases to ~ 6 MeV when 
some nonzero modes quantum corrections are included [9] . Some states with positive parity 
and unusual connection between isospin and angular momentum are predicted also, which 
can reveal themselves as supernarrow radiatively decaying dibaryon [10], as well as dibaryons 
with negative parity. The experimental indications on the existence of such states have 
been obtained recently [11, 12]. Another issue is the binding energy of strange hypernuclei 
with atomic number up to ~ 16 which can be described qualitatively under some natural 
assumptions, and predictions of more bound charmed or beautiful hypernuclei [13]. 

Here we shall discuss the nuclides with baryon number A = B < 32: mass splittings of 
isotopes with different N, Z, including neutron rich nuclides with the difference N — Z up to 
10 and 11, N, Z are the numbers of neutrons and protons in the nucleus, N + Z = A. Some 
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extrapolation and modification of the approach applied previously [4, 5, 13] allows to make 
certain predictions and conclusions also in this case. In particular, the symmetry energy term 
in the Weizsacker — Bethe — Bacher mass formula for nuclei [14], E sym = b sym (N — Z) 2 /2A 
can be described successfully (b sym ~ 50 MeV). The presence of this term (more correctly, 
it should be called the asymmetry energy since it appears due to difference of iV and Z 
numbers) combined with the energy of Coulomb repulsion of protons, Eq ~ 0.6Z 2 e 2 / Rz, 
leads to the observed excess of neutrons in nuclei (see, e.g. [15, 16]). Within the chiral soliton 
approach this symmetry energy originates from isotopical rotations of multiskyrmions and 
has the form E sym = 1(1 + l)/(20/), with isospin / = (N — Z)/2 for the states of lowest 
energy (we shall consider mainly the case N > Z, as it follows from presence of the Coulomb 
repulsion). Within the conventional approach this term appears due to dependence of the 
Fermi-motion kinetic energies ez, e^v of protons and neutrons on N and Z numbers [15]: 

2 

Osym — -{tZ + £N), 

ez — kp(p) 2 /2M, with Fermi-momentum of protons in the nucleus kp(p) = (3ir 2 Z /Vp) 1 / 3 , 
V p being the volume occupied by protons in the nucleus, and similarly for ejy (V n can be 
different from V p ). Numerically this gives for heavy nuclei b sym ~ 24 MeV if the radius 
of the nucleus (for protons and neutrons) Ra = 1.12 A 1 / 3 fm, and approximately the same 
contribution give potential terms [15]. The presence of the linear in isospin / term shifts 
slightly the balance of energies E syra and E c towards larger values of N, in better agreement 
with data, especially for smaller A. 

The question, if the neutron rich light nuclides may exist, which we address here 
as well, seems to be more viable now when the beams of nuclides become available, and 
the possibilities for production and observation of such states have been increased [17—20]. 
After short theoretical introduction in the next section we discuss qualitatively situation 
with the mass splittings of relatively light nuclear isotopes (Section 3), give the results for 
differences of binding energies of isotopes with integer isospin (Section 4) and half-integer 
isospin, odd atomic numbers (Section 5). Some predictions for binding energies of neutron 
rich nuclides are given in Sections 4 and 5. Discussion of our method and results is presented 
in concluding Section 6. 



2. Theoretical background 

The brief description of the chiral soliton approach and some basic formulas are 
necessary. The SU(2) skyrmions are described by the Lagrangian density of the effective 
field theory, depending on the SU(2) matrix U = f + i(ri/i + T2/2 + ^"3/3)- The Lagrangian 
in its "minimal" form, suggested by Skyrme [1], is 

F 2 1 F 2 m 2 

L = TrL^ + ^TrfVM 2 + -^Tr(U + tf-2) (1) 

with L M = d,JjW = Hp ■ r. 

The sixth order in chiral derivatives term also can be included into consideration 

[21] i): 

U = -c 6 Tr([L M ^][L^][L 7 i/]) 

1 ) Probably, one of the first attempts to include the sixth order term was made in [7] where the bound B = 2 
torus-like configuration was found, similar to the case of fourth order, or Skyrme term. 



with dimensional constant cq. The important property of the Lagrangian (1) is that it is 
proportional to the number of colors iV c of the underlying quantum chromodynamics - 
theory of colored quarks and gluons and their interactions [22]. The effective field theory 
described by Lagrangian density (f ) is being built from this QCD. 

The baryon, or winding, number is the fourth component of the Noether current 
generated by the Wess-Zumino term in the action [22], and in these notations equals to 



B 



^- J Tr^L^tfV = - JL / (l^ydV, (2) 



(abc) denotes the mixed product of vectors a, b, and c. We shall rewrite it in less conven- 
tional but transparent form 

B = ~ J 000f 3 )/f o dPr = 

= ~hl 6{f2 ° + f " + f * + f " - l )^ B ^ B ^W r - ( 3 ) 

Three functions are introduced usually to describe ff f = cp, fk = spUk, where the 
components of the unit vector n are rt\ = s Q cg, n 2 = s^s^, n 3 = c Q , cp = cosF, sp = sinF, 
s a = sin a, etc. Then (df 1 df 2 df 3 )/f = s 2 F s a (dFdad/3), and the winding number 

B = / s 2 F s a I[(F, a, P)/(x, y, z)}d 3 r (4) 

with I[(F, a, /3)/(x, y, z)\ - Jacobian of the transformation from variables r = x, y,z to 
F, a, (3. Since the element of the unit 3-dimensional sphere dS 3 = s 2 F s a dF dadfi, and 2-k 2 is 
the surface of the unit sphere, Eq. (4) shows how many times the sphere S 3 (homeomorphic 
to SU(2)) is covered when integration over d 3 r is made, so it is the degree of the map 
R 3 -> SU(2). 

The static energy of arbitrary SU (2) skyrmion in these notations is 
M = J |^[(9>) 2 + s 2 F {0af + slim) 2 )} + 

+ ^([dFx da] 2 + s 2 a [dF x d(3} 2 + s 2 F s 2 a [da x <9/3] 2 ) + m 2 F 2 (l - c F )/A + 

+ 9Qc 6 s F s 2 a 0Fdad(3) 2 ^d\ (5) 

[a x b] means the vector product of a and b. The mass (5) is proportional to the number 
of colors of underlying QCD, i.e. it is large at large N c . 

The rational map approximation (RM) proposed in [23] is very effective for descrip- 
tion of multibaryons with large enough baryon numbers. The direct numerical calculation of 
configurations of minimal energy [24] provides the results very close to those obtained within 
the RM approximation. The profile function F within the RM approximation depends on 
variable r only, and the unit vector n — on angular variables 9, 0: n\ = 2Rei?/(l + |-R| 2 ), 
n 2 = 2Imi?/(l + \R\ 2 ) and n 3 = (1 - |i?| 2 )/(l + \R\ 2 ), R(6, 0) being a rational function of 
variable z = tan(0/2) exp(i0) defining the map of degree M of the 2-dimensional spheres 
S 2 — > S* 2 . We shall use the following notations, introduced in [23]: 
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dfl = smOdOdcf) = 2idzdz/{l + \z\ 2 ) 2 . For configurations of lowest energy Af = B and 
F(0) — F(oo) = 7T. The classical mass of the configuration for the RM ansatz can be written 
as [23, 21] 
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JL J | r '2 F /2 + 2Bs 2 p[1 + (1 _ A)F / 2] +1^1 + \(F' 2 - 1)] }rfr + M.t., (8) 



where the distance r is measured in units 2/(F n e'), e' = e^l — A, A defines the weight of 
the sixth order term in the Lagrangian according to the relation A/ (1 — A) 2 = 48c6-F 2 e 4 . If 
A = 0, c 6 = and we obtain original variant of the Skyrme model (SK4 variant), A = 1 
corresponds to the pure SK6 variant. In this case relation takes place: e' = l/(48c 6 F 2 ) 1 / 4 . 
M.t. denotes the contribution of the mass term. 

The energy of the system depending on the angular velocities of rotations in SU(2) 
isospin collective coordinate space can be written in such form: 
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The angular velocities of rotation in the isospace are defined in standard way: A^A = 
The functions Ui are connected with the body fixed angular velocities of SU (2) rotations by 
means of transformation: 

(V ■ T = U ] U) tU -u t, (10) 

withl^ + l^ + l 2 = ( ( 9/ ) 2 + (9/ 1 ) 2 + (9/ 2 ) 2 + (9/ 4 ) 2 and (I; 2 =4[w 2 f 2 -(wf) 2 ]. The coefficients 
in the quadratic form (9) in oJiUOk define the SU(2) (isotopic) moments of inertia of arbitrary 
systems of SU(2) skyrmions. In the axially symmetrical case and in RM approximation 
they are not complicated and were written explicitly in [8, 25, 26]. 

Some further complication is necessary because one should introduce another set 
of collective coordinates [8] to take into account the usual space rotations which we shall 
describe by angular velocities i = 1, 2, 3. The part of the Lagrangian describing both 
kinds of rotations can be written then [8]: 

1 1 

L rot = -UijUiUj + -Vijfliflj - WijUiQj, (11) 

where tensors Uy, V^, and are defined by the configuration of lowest energy for each 
value of B, i.e. they are functions of F, a, and (3. 

The body fixed isospin and spin angular momentum operators K\ = dL rot /duii and 
L\ = dL Iot /dQi are linear combinations of u>i and Three diagonal moments of inertia and 
three off-diagonal define each of three rotation terms in (11) in most general case. In the 
case of axially symmetric systems we obtained four different diagonal moments of inertia: 
@i = @2 = @at; ©3? @_r and 6 int [8] 2 \ The static masses of solitons and momenta of 
inertia were calculated previously for B < 32 starting from rational map ansatz [23, 24] 
and using the variational minimization program [25, 26], see Table 1. It should be noted 
that expressions (1) — (9) contain all information which is necessary for the description of 
arbitrary SU(2) skyrmions, in such minimal form. It was sufficient for the fitting of the 

2 ^For the hedgehog- type configurations the isospin and usual space tensors of inertia coincide since iso- 
and usual space rotations are identical for the hedgehogs, and the configuration is described by one common 
moment of inertia. 



mass differences of the nucleon and A isobar [2, 3] and for the description of many basic 
properties of light nuclei [4, 5]. 

The moments of inertia, isotopical and orbital, necessary for our purpose here, are 
given by the expressions: 
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= T / + ^ [ (1 " A) + X i ) + 4AJ/ 

The inequalities take place for any value of A [25, 26]: 
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which ensures that the quantum correction due to usual space collective rotations is always 
smaller (even much smaller at large A = B) than isotopical quantum correction. The 
moments of inertia (12), (13) are proportional to the number of colors N c , similar to the 
classical mass of solitons, therefore, the rotational quantum corrections are proportional to 
~ 1/iVc, i.e. they are small, and the whole consideration becomes selfconsistent at large N c . 
The fact that N c = 3 in our world makes the real life somewhat more complicated. 

Considerable simplification takes place because of the symmetry properties of many 
multiskyrmion configurations, beginning with 5 = 3,7 and some other configurations found 
in [23, 24]. It turned out that all three tensors of inertia in (11) at large B approximately 
are proportional to the unit matrix, i.e. £4, ~ uSij, Vij ~ v5 i j 1 Wij ~ wSij. Moreover, in 
many cases of interest the interference tensor of inertia is negligibly small, w <^iu and ti)<« 
[26, 27]; the relative magnitude of this interference inertia decreases with increasing baryon 
number, and we shall neglect it for numerical estimates. 

The mass splitting between nucleon and A (1232) — usually one of the first quantities 
fitted within the chiral soliton approach — is described almost equally for different choices 
of the parameters used in the literature [2, 3, 26] (it scales like ~ F^e 3 approximately), 
but dimensions of solitons are different since they scale like l/(F n e). One of the possible 
choices of the parameters of the model is [28], F n = 186 MeV, e = 4.12, and we shall 
use these numbers for the "nucleon" variant of the model. Dimensions of multiskyrmions, 
the classical field configurations, are much smaller than dimensions of observable quantized 
states, i.e. nuclei. The latter are reproduced only when the vibration and breathing motions 
in the solitons are taken into account, as it was shown recently for the deuteron in [9]. The 
numerical values of inertia are given in Table 1 for the SK4 (e = 4.12 and e = 3.0, marked 
with *) and SK6 variants of the model, e' = 4.11 and 2.84, also marked with *. The values 
e = 3.0 and e' = 2.84 for the SK6 variant allow to describe well the isotopical splittings of 
nuclei with atomic numbers between ~ 10 and 32, see discussion in the following sections. 
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Table 1. The moments of inertia of multiskyrmions in the SK4 variant of the model e = 4.12 
and 3 (6j(SK4)*), and for the SK6 variant of the model, e' = 4.11 and 2.84 (6j(SK6)*), in 
GeV" 1 



3. Multiplets of nuclear isotopes 

To get an idea how the chiral soliton approach can describe the existing data on 
the mass splittings of the known nuclear isotopes, we consider the mass differences of some 
known nuclei with atomic number beginning with A = 8. 

Consider first the mass difference between isosinglet 8 Be, J p = + , binding energy 
(b.e.) e = 56.5 MeV, and the components of the isotriplet 8 Li- 8 B, J p = 2 + , b.e. 41.28 MeV 
and 37.74 MeV (the experimental data on total binding energies of nuclei are taken mainly 



from paper [29]). Using the above formula, we obtain easily 



AM(A = 8)0,1 = [M( 8 Li)+M( 8 B)]/2-M( 8 Be) = -— + -_. (15) 

In the average of masses of isotopes 8 Li and 8 B the number of protons and neutrons is the 
same as in 8 Be, and the average Coulomb energy is also the same with quite good accuracy. 
The theoretical value ~ 40 MeV for the SK4 variant of the model should be compared with 
17 MeV from the data. Only collective motion is taken into account in (15), and peculiarities 
of nucleon-nucleon interaction are totally neglected, e.g. the phenomenologically introduced 
so called pairing energy [15, 16] which increases the binding by ~ 4 MeV for 8 Be and 
decreases by same value for 8 Li and 8 B. As a result we should compare 40 MeV and 
~ 10 MeV. So, in this case we can speak only about qualitative agreement. 

For A = 12 there is isosinglet nucleus 12 C with J p = + and b.e. 92.16 MeV and 
isotriplet components 12 B and 12 N with J p = 1 + , b.e. 79.58 and 74.04 MeV, 

AM(A = 12)o,i = ^~ + (16) 

W,12 W,12 

which is 22.2 MeV to be compared with 15.35 MeV experimentally. The pairing energy in 
this case also makes the disagreement even worth by 5—6 MeV. 

For A = 16 isosinglet and scalar nucleus 16 0, b.e. 127.62 MeV should be compared 
with nucleus 16 N, J p = 2", b.e. 117.98 MeV and 16 F, b.e. 111.41 MeV, if we assume that 
the last one also has J p = 2~ . In this case 

AM{A = 16)o,i = 7^ + 7^-, (17) 

numerically it is 18.55 MeV to be compared with 14.69 MeV, or with ~ 10 MeV when 
pairing energy is added and subtracted in proper way. 

For the case of A = 10 we have isosinglet 10 B, J p = 3 + , e = 64.75 MeV and compo- 
nents of isotriplet 10 Be and 10 C, J p = 0+, b.e. 64.98 and 60.32 MeV. Theoretical difference 
of masses is 

AM{A = 10)o,i = ^-~ (18) 

which is about 11 MeV, to be compared with experimental value 2.1 MeV. The pairing 
energy of the order of 3—4 MeV decreases the binding of 10 B (N = Z = 5) and increases the 
binding of 10 Be and 10 C, so we should compare 11 MeV and 8—10 MeV. 

To summarize, we can state that the agreement between data and theory is only 
qualitative for small A — B, if we take the value of the model parameter e = 4.12 which 
allowed to describe the mass difference between nucleons and A(1232) isobar (we shall call 
this as "nucleon" variant of the model), which improves however with increasing atomic 
number, as can be seen for A = 12, 16. 

It is clear that the "nucleon" variant of the model cannot be good in the case of nuclei: 
the chiral field configurations given by rational map ansatz, although confirmed by direct 
numerical computation, represent the "sceleton" configurations which differ from realistic 
distribution of nucleons inside of nuclei. The nonzero mode motion should be taken into 
account — vibration, breathing, etc. to get more realistic distribution of matter inside of 
multibaryons. Technically, this problem is very complicated. Up to now, the mathematical 
analysis of these modes was made for baryonic numbers up to A = 7 [30], and some numerical 
estimates were made for the deuteron (A = 2) [9], as it was mentioned in Introduction. For 
each value of A we can improve the agreement of theoretical estimates with data if we allow 



the only parameter of the model, Skyrme constant e, to vary — to decrease, in fact, in 
order to make the dimensions of multiskyrmion and moments of inertia greater. Effectively, 
such approach allows to take into account the vibration and breathing of separate parts of 
multiskyrmions. Evidently, for each value of atomic number we can get perfect agreement 
if e = e(A). It is, however, more instructive to describe a number of nuclei with the same 
value of constant e. The optimal value of the constant is e = 3.0 (SK4 variant) and e = 2.84, 
as analysis of data shows (next section) . 

With rescaled value of the constant e we obtain (see Table 1): AM (A = 8) = 
16.8 MeV; AM(A = 12) = 9.5 MeV and AM(A = 16) = 7.9 MeV, and for A = 12 and 16 
there is now agreement with data within 2—3 MeV. As we shall see below, quite satisfactory 
agreement with data can be obtained when appropriate choice of nuclear isotopes is made, 
for A between ~ 10 and ~ 30. 



4. Even baryon numbers, integer isospins 

It is possible to obtain more adequate comparison of the model predictions with data, 
comparing the differences of masses of nuclei with even isospins, odd isospins, and similar 
differences for half-integer isospins (next section). As it was mentioned in previous section, 
the nucleons pairing energy should be added to the binding energy, which equals to A when 
both N, Z are even, when A is odd, and —A when both N, Z are odd (see, e.g. [15]). 
This pairing energy decreases with increasing atomic number; it means that peculiarities 
of nucleon interaction become less important for large A. In the differences of masses of 
isotopes with the same pairing energy this specific contribution is cancelled, and one can 
hope that the collective motion effects — which are mainly taken into account within the 
chiral soliton approach — play the dominant role. 

Formulas for energies differences can be correctly applied to the same nucleus' quan- 
tum states with different isospin. However, experimental data for excited states with higher 
isospin are lacking often. Nevertheless, using the isotopical invariance we can obtain binding 
energies of this level subtracting the Coulomb energy difference from binding energies of the 
components of isotopical multiplet with neutron excess. 

In Tables 2—10 the Coulomb energies are calculated according to [16] (finite-range 
liquid-drop model): 

3 e 2 Z 2 3 / 3 \ 2 / 3 e 2 Z 4 / 3 1 /145\ r 2 p e 2 Z 2 

C ~hr A^ * 4I2J roA 1 ' 3 8 I 48 ) r§A ' { } 

where nuclear-radius constant r = 1.16 fm, proton root-mean-square radius r p = 0.80 fm. 
The first term in (19) is Coulomb energy for a homogeneously charged, diffuse-surface nu- 
cleus to all orders in the diffuseness constant a dcn , the second is Coulomb exchange correc- 
tion, and the third is proton formfactor correction to the Coulomb energy. The constant 
B 3 = B 3 (tqA 1 / 3 /a dcn ) , is normalized so that for zero diffuseness a den (range of Yukawa 
function exp(— r/a^)) B 3 — 1. For spherically symmetrical nuclei it can be calculated an- 
alytically and presented as a sum £> 3 = 1 — 5y 2 + here y = ad cn /(^o^ 1 ^ 3 )? more details 
can be found in [16]. 

The result provided by (19) for intermediate values of A (practically, between A ~ 25 
and A ~ 50) can be described by simple formula 



(20) 



where Q is nucleus's charge, tq = 1.48 fm [14]. The effective Coulomb radius in (20) is greater 
than r c given, e.g. in [15], r c = 1.24 fm. The formula (20) is valid for uniformly charged 
sphere of radius R = tqA 1 ^ and total charge Q. For other types of charge distribution 
inside of nuclei the above formula for Eq is modified slightly. For example, for thin shell- 
like distribution suggested by rational map approximation for multiskyrmions, the Coulomb 
energy equals to 

e - - \S*< (21) 

which is smaller than in the case of uniformly charged sphere by 17% only. 

For the nuclei considered in Tables 2—5 and 7—9, the Coulomb energy increases from 
zero for nucleus 6 H to 43.9 MeV for nucleus 32 S. 



A 


cxp 


cxp 


Ae c 


A cxp 
^ e 02 


A th 
^ Ae 02 


A th* 
^ e 02 


6 Li- 6 H 


32.0 


5.8 


1.1 


27.3 


112.9 


47.9 


8 Be- 8 He 


56.5 


31.4 


2.0 


27.1 


88.5 


37.5 


io B _io Li 


64.8 


45.3 


2.9 


22.4 


72.5 


30.8 


i2 C _i2 Be 


92.2 


68.7 


3.8 


27.3 


61.8 


26.3 


14 N _14g 


104.7 


85.4 


4.6 


23.9 


53.5 


22.8 


16 _16 C 


127.6 


110.8 


5.4 


22.2 


47.5 


20.3 


18p_18 N 


137.4 


126.5 


6.1 


17.0 


42.7 


18.2 


20 Ne- 2U O 


160.6 


151.4 


6.8 


16.0 


38.7 


16.6 


22 Na- 22 F 


174.1 


167.7 


7.5 


13.9 


35.6 


15.2 


24 Mg- 24 Ne 


198.3 


191.8 


8.2 


14.7 


32.9 


14.1 


26 Al- 26 Na 


211.9 


208.2 


8.9 


12.6 


30.5 


13.1 


28 Si- 28 Mg 


236.5 


231.6 


9.5 


14.4 


28.5 


12.3 


30p_30 i y 


250.6 


247.8 


10.2 


13.0 


26.8 


11.5 


32 S- 32 Si 


271.8 


271.4 


10.8 


11.2 


25.3 


10.9 



Table 2. The binding energies differences for isotopes with isospin 1 = and 2 for the orig- 
inal variant, e = 4.12, and for the variant with rescaled constant, e = 3 (numbers with the *) 



The pairing energy has different sign for nuclei with A = 6, 10, 14,... and A = 
8, 12, 16, but anyway, it is cancelled in the binding energies differences we consider here. 

To obtain the binding energy of isotopes with isospin I = we take the known 
binding energy of ground states of nuclei with N = Z and add the energy of the Coulomb 
repulsion given by (19), (20). In the case of / = 2 we take the ground state of nucleus with 
N = Z + 4 which is supposed to be the component of / = 2 isomultiplet with the isospin 
projection J 3 = —2, and add the Coulomb repulsion energy as well. The quantity Ae exp 
shown in all Tables includes the Coulomb correction. 

The energy of usual space rotation of nucleus as a whole, equal to J(J + l)/20j, 
is small for large enough nuclei since Oj is large (see Table 1), and is omitted within this 
approach. Another reason to omit this correction is that the value of J is not known in 
some cases. 

It follows from Table 2 that for light nuclei, A between 6 and 10, there is only 
qualitative agreement between data and theoretical results, even with rescaled constant e. 
For atomic numbers between 12 and 32 the agreement is satisfactory and even good in some 
cases (see column with Ae th> *. By this reason we shall not consider further small values of 
A. 



A 


cxp 
e l 


cxp 
e 3 


Ae c 




A th 
AAfc 13 


A th* 
AAfc 13 


10 Be- lu He 


65.0 


30.3 


2.0 


36.7 


120.8 


51.3 


i2 B- i2 Li 


79.6 


44.4 


2.9 


38.1 


103.0 


43.8 


14 C- 14 Be 


105.3 


70.0 


3.7 


39.0 


89.1 


38.0 


16 N- 16 B 


118.0 


88.2 


4.5 


34.3 


79.1 


33.8 


18 0- 18 C 


139.8 


115.7 


5.2 


29.3 


71.2 


30.4 


20p_2U N 


154.4 


134.2 


6.0 


26.2 


64.5 


27.6 


22 Ne- 22 


177.8 


162.0 


6.7 


22.5 


59.3 


25.4 


24 Na- 24 F 


193.5 


179.1 


7.4 


21.8 


54.8 


23.5 


26 Mg- 26 Ne 


216.7 


201.6 


8.1 


23.2 


50.9 


21.9 


28 Al- 28 Na 


232.7 


218.4 


8.7 


23.0 


47.5 


20.4 


30 Si _30 Mg 


255.6 


241.6 


9.4 


23.4 


44.6 


19.2 


32p_32^] 


270.9 


259.2 


10.0 


21.7 


42.2 


18.2 



Table 3. The binding energies differences for isotopes with isospin 1 = 1 and 3 for the 
original variant, e = 4.12, and for the variant with rescaled constant, e = 3 

In the case presented in Table 3 we take the binding energies of isotopes with iV = 
Z + 2 and iV = Z + 6 (ground states) which are considered as components of isomultiplets 
1 = 1, h = — 1 and 1 = 3, Is = —3, Coulomb energy is included into consideration as well. 
The agreement of theoretical estimates with data is quite satisfactory for A between 14 and 
32, qualitative agreement takes place for A = 10, 12 as well. 



A 


cxp 

e o 


cxp 
e 4 


Ae c 


A cxp 
AAt 04 


A th 
AAfc 04 


A th* 
AAfc 04 


16 Q _16 Be 


127.6 




9.0 




158.3 


67.5 


i«F- 18 B 


137.4 


89.0* 


10.5 


59.4 


142.3 


60.8 


2U Ne _2U C 


160.6 


119.2 


12.0 


53.4 


129.1 


55.2 


22 Na- 22 N 


174.1 


140.0 


13.4 


47.5 


118.6 


50.8 


24 Mg- 24 


198.3 


168.5 


14.8 


44.6 


109.5 


47.0 


26 A1- 26 F 


211.9 


184.5 


16.2 


43.6 


101.8 


43.7 


28 Si- 28 Ne 


236.5 


206.9 


17.5 


47.1 


95.1 


40.9 


3Up_3U Na 


250.6 


224.9 


18.8 


44.5 


89.3 


38.4 


32 S- 32 Mg 


271.8 


249.7 


20.1 


42.2 


84.4 


36.3 



Table 4. The binding energies differences for isotopes with isospin 1 = and 4 for the 
original variant, e = 4.12, and for the variant with rescaled constant, e = 3 

For I = A the ground states of isotopes with N = Z + 8 are taken as I 3 = —4 
components of corresponding isomultiplets. The agreement of data and theoretical values 
for atomic numbers between 18 and 30 is good or satisfactory, according to Table 4. For the 
nucleus 16 Be we obtain the b.e. about 60 MeV, so it is unstable: it can decay into isotope 
14 Be (b.e. about 70 MeV) and 2 neutrons. This result can be considered as prediction of our 
approach. For the 18 B isotope our prediction for the b.e. is 87.1 MeV 3 \ in good agreement 
with extrapolation result 89 MeV of [29], marked in Table 4 with t. 

Results presented in Tables 2—4 are illustrated in Fig. 1. 

3 )The value 87.1 MeV = (137.4-60.8 + 10.5) MeV, where 60.8 MeV is the difference of quantum corrections 
to the multiskyrmion, rescaled variant, and 10.5 MeV is difference of Coulomb energies 



Ae 70 




1 1 1 1 1 1 1 1 1 1 1 ' 1 1 1 

6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 A 

Fig. 1. The binding energies differences (in MeV) for isotopes with even atomic numbers 
for the SK4 variant with rescaled constant, e = 3.0 (black points connected with solid lines 
— experimental data, dashed lines — model calculations) 



A 


exp 
e 


exp 
e 2 


Ae c 


A exp 
^ e 02 


A SK6 
AAfc 02 


A SK6,* 


8 Be- 8 He 


56.5 


31.4 


2.0 


27.1 


110.1 


38.8 


iu B _iu Li 


64.8 


45.3 


2.9 


22.4 


91.0 


32.1 


12 C- 12 Bc 


92.2 


68.7 


3.8 


27.3 


78.0 


27.5 


14 N _14 B 


104.7 


85.4 


4.6 


23.9 


67.7 


23.9 


16 0- 16 C 


127.6 


110.8 


5.4 


22.2 


60.3 


21.3 


i«F- 18 N 


137.4 


126.5 


6.1 


17.0 


54.3 


19.2 


2U Ne _20 Q 


160.6 


151.4 


6.8 


16.0 


49.3 


17.5 


22 Na- 22 F 


174.1 


167.7 


7.5 


13.9 


45.3 


16.1 


24 Mg- 24 Ne 


198.3 


191.8 


8.2 


14.7 


41.9 


14.9 


26 Al- 26 Na 


211.9 


208.2 


8.9 


12.6 


38.9 


13.8 


28 Si- 28 Mg 


236.5 


231.6 


9.5 


14.4 


36.4 


13.0 


30p_30 i y 


250.6 


247.8 


10.2 


13.0 


34.1 


12.2 


:!2 S :i2 Si 


271.8 


271.4 


10.8 


11.2 


32.3 


11.5 



Table 5. The binding energies differences for isotopes with isospin 1 = and 2 for the SK6 
variant, e' = 4.11, and for the SK6 variant with rescaled constant, e' = 2.84 

In Table 5 we present the comparison of data for differences of binding energies of 
isotopes with isospin and 2 (as in Table 2) and results of the SK6 variant of the model, for 
two values of effective constant e. The rescaling of the constant e' to the value 2.84 allows 
to make agreement of theoretical prediction and data quite satisfactory, beginning with 
atomic number A = 12, similar to the SK4 variant of the model. We conclude therefore that 
there is no difference of principle between both variants of the model, and will make further 



comparison with the results provided by the SK4 variant of the model. Our predictions for 
the binding energies of nuclides with highest values of isospin (Tables 6, 10) will be made 
for both variants of the model. 

In view of successful description of binding energies of many isotopes, it is possible 
to make predictions of binding energies for neutron rich nuclides which have greater isospin 
than we have considered before, e.g. / = 5 or N = Z + 10. 



A 


14/ @SK4 


Ae c 


C SK4 
e /=5 


C SK6 
e /=5 




cxp 

e/=5 


l8 Be 


85.4 


8.8 


63 


59 






2U B 


77.3 


10.3 


87 


83 






22 C 


71.1 


11.8 


118 


114 






24 N 


65.7 


13.2 


141 


138 






26Q 


61.1 


14.6 


170 


167 


168 




28p 


57.1 


15.9 


191 


188 


186 




30 Ne 


53.8 


17.2 


219 


216 




212 


32 Na 


50.9 


18.5 


238 


235 




231 



Table 6. Predictions for the binding energies of neutron rich nuclides with N = Z + 10 or 
I = 5 for the SK4 and SK6 variants of the model. The values are the results of extrapo- 
lation of [29] 

The binding energy of neutron rich nuclides with N = Z + 10 shown in Table 6 is 
calculated according to the following formula: 

14 

£n=z+io — £n=z+2 — V Aec, (22) 

<Oa=N+Z 

where e^=z+2 is the experimentally known binding energy of the isotopes with isospin I — 1, 
the number 14 is half of difference of Casimir operator /(/ + 1) for 1 = 5 and 1=1, Aec is 
the difference of the repulsive Coulomb energies for nuclides with the same A = N + Z and 
different Z, calculated with formula (19). Predictions of two variants of the model do not 
differ much, the maximal difference is about 4 MeV for lighter nuclides. 

5. Odd baryon numbers, half-integer isospin 

The consideration similar to that made for integer isospin in previous section can 
be performed for half-integer isospins, i.e. for nuclei with odd atomic numbers. The phe- 
nomenologically introduced pairing energy is absent in this case. 



A 


cxp 
C l/2 


cxp 
fc 5/2 


Ae c 


A exp 
^1/2,5/2 


A th 

^1/2,5/2 


A th* 
^1/2,5/2 


9 Be- 9 He 


58.2 


30.3 


2.0 


29.9 


105.8 


45.0 


n B- n Li 


76.2 


45.6 


2.9 


33.5 


88.6 


37.7 


13 C- 13 Be 


97.1 


68.1 


3.7 


32.7 


76.8 


32.7 


15 N- 15 B 


115.5 


88.2 


4.5 


31.8 


66.9 


28.6 


17 _17 C 


131.8 


111.5 


5.3 


25.6 


60.4 


25.8 




147.8 


132.0 


6.0 


21.8 


54.1 


23.2 


21 Ne- 21 


167.4 


155.2 


6.8 


19.0 


49.4 


21.2 


23 Na- 23 F 


186.6 


175.3 


7.5 


18.8 


45.5 


19.5 


25 Mg _25 Ne 


205.6 


196.0 


8.2 


17.8 


42.3 


18.1 


27 Al- 27 Na 


225.0 


214.9 


8.8 


18.9 


39.4 


16.9 


29 Si- 29 Mg 


245.0 


235.3 


9.5 


19.2 


36.9 


15.9 


31p_31 A1 


262.9 


255.0 


10.1 


18.0 


34.7 


14.9 



Table 7. The binding energies differences for isotopes with isospin 7 = 1/2 and 5/2 for the 
original variant, e = 4.12, and for the variant with rescaled constant, e = 3 

In Table 7 we take the binding energies of nuclei with N = Z + 1 for isospin 1 = 1/2 
and N = Z + 5 for / = 5/2. For the rescaled variant of the model the quantitative agreement 
with data is quite good for atomic numbers between 13 and 25, qualitative agreement takes 
place for all values of A presented in Table 7. 



A 


cxp 
C 3/2 


cxp 
C 7/2 


Ae c 


A cxp 
AAt 3/2,7/2 


a th 

^ e 3/2,7/2 


A th* 
AAt 3/2,7/2 


15 C- 15 Be 


106.5 




3.7 




100.4 


42.8 


17 N- 17 B 


123.9 


89.6 


4.4 


38.7 


90.6 


38.7 


19 _19 C 


143.8 


115.8 


5.2 


33.2 


81.2 


34.7 


2ip- 21 N 


162.5 


138.8 


5.9 


29.6 


74.1 


31.8 


23 Ne _23Q 


183.0 


164.8 


6.6 


24.8 


68.2 


29.2 


25 Na- 25 F 


202.5 


183.5 


7.3 


25.3 


63.4 


27.2 


27 Mg _27 Ne 


223.1 


203.0 


8.0 


28.1 


59.0 


25.4 


29 Al- 29 Na 


242.1 


222.8 


8.7 


28.0 


55.4 


23.8 


31 Si- 31 Mg 


262.2 


244.0 


9.3 


27.5 


52.1 


22.4 



Table 8. The binding energies differences for isotopes with isospin / = 3/2 and 7/2 for the 
original variant, e = 4.12, and for the variant with rescaled constant, e = 3 

The nuclei with N = Z + 3 and N = Z + 7 in Table 8 correspond to the states 
with isospin / = 3/2, 1% = —3/2 and / = 7/2, 1% = —7/2. Good agreement between data 
and theory (rescaled variant) takes place between A = 17 and 27. The nucleus 15 Be is not 
observed yet, so, the b.e. of this nucleus 63.7 MeV is prediction of the model. It is not 
stable since it can decay into 14 Be (b.e. 70 MeV) plus neutron. 

Satisfactory agreement between available data and theoretical result (rescaled variant 
of the model with Skyrme constant e = 3.0) takes place in the case of isotopes with isospin 
1/2 and 9/2 as well, see Table 9. The numbers with t in the 3-d and 5-th column are results 
of extrapolation made in [29] . 

The nucleus 17 Be, similar to 15 Be, has not been observed yet, so, the b.e. 63.4 MeV 
which follows from Table 9 also is prediction of the chiral soliton approach. This isotope is 
not stable since the b.e. of 14 Be, 70 MeV, is greater. For the isotope 19 B the prediction of 
our approach is 88.7 MeV, in good agreement with the extrapolation result 90.1 MeV shown 
in Table 9. Satisfactory agreement between our prediction and extrapolation of paper [29] 
takes place for isotopes 21 C, 23 N, and 25 0. 



A 


cxp 
C l/2 


cxp 
e 9/2 


Ae c 


A cxp 
^1/2,9/2 


A th 

ZAC l/2,9/2 


A th* 
ZAC l/2,9/2 


17 0- 17 Be 


131.8 




8.9 




181.3 


77.3 


iyp- iy B 


147.8 


90.1* 


10.4 


68. 1 1 


162.4 


69.5 


21 Ne- 21 C 


167.4 


118.8 7 


11.9 


60.5 T 


148.3 


63.5 


23 Na _23 N 


186.6 


142.2* 


13.3 


57.7* 


136.4 


58.5 


25 Mg- 25 


205.6 


168.4* 


14.7 


51.9* 


126.8 


54.4 


27 A1 _27p 


224.9 


185.8 


16.0 


55.1 


118.1 


50.7 


29 Si- 29 Ne 


245.0 


208.2 


17.4 


54.2 


110.7 


47.6 


3i P- 3i Na 


262.9 


228.9 


18.7 


52.7 


104.2 


44.8 



Table 9. The binding energies differences for isotopes with isospin 1 = 1/2 (i.e. N = Z +1) 
and 9/2 (N = Z + 9) for the original variant, e = 4.12, and for the variant with rescaled 




Fig. 2. The binding energies differences (in MeV) for isotopes with odd atomic numbers 
for the rescaled SK4 variant of the model, e = 3.0 (black points connected with solid lines 
— experimental data, dashed lines — model calculations) 



constant, e = 3.0 



The results for odd atomic numbers (half integer isospins, Tables 7—9) are presented 
also in Fig. 2. Similarity to behavior of data and calculation results shown in Fig. 1 is 
evident. 

In Table 10 we present also predictions of the chiral soliton approach for nuclides 
with TV = Z + 11, or isospin /= 11/2 for the SK4 and SK6 variants of the model. 



A 


16/ @SK4 


Ae c 


,SK4 
e ll/2 


C SK6 
e ll/2 


J 

C ll/2 


17 Li 


103.2 


7.2 


28 


23 




19 Be 


92.5 


8.7 


60 


56 




2l B 


84.7 


10.2 


88 


83 




n C 


78.0 


11.7 


117 


113 




25 N 


72.4 


13.1 


143 


139 




27Q 


67.5 


14.4 


170 


166 




29p 


63.5 


15.8 


194 


191 


187 


31 Ne 


59.7 


17.1 


219 


216 


212 



Table 10. Predictions for the binding energies of nuclides with / = 11/2 or N = Z + 11 
calculated from nuclides with 1 = 3/2 (N = Z + 3) 

For the nuclides 29 F and 31 Ne we obtain fair agreement with the extrapolation result 
obtained in [29], and the SK6 variant prediction is more close to it than the SK4 one. 
Experimental data for this case and A < 32 are still lacking. 

6. Conclusions 

The estimates of the mass splittings of nuclides with baryon numbers up to 32 are 
made using the properties of the bound states of skyrmions obtained in the rational map 
approximation. The results are in impressive agreement with data, provided the rescaling 
of the Skyrme constant e is made — decrease by about 25—30% — which effectively takes 
into account the role of nonzero modes, vibration and breathing. This nonzero modes 
motion brings the rational map skyrmions, classical configurations of the one-shell type, in 
better agreement with observed nuclear matter distribution: it increases the dimensions of 
configuration and, therefore, the moments of inertia. Without rescaling of the constant e 
the mass splittings of light nuclear isotopes are overestimated more than twice. It is known 
that calculations of nonzero modes contribution to the energy (mass) of multiskyrmion are 
closely related to calculations of the so called Casimir energy and/or loop corrections to the 
soliton mass [31] which are of the order of N® ~ 1. The estimate of this contribution was 
made for the B = 1 hedgehog-type skyrmion [31]. Anyway, these unknown contributions 
cancel out in the differences of energies of quantized states we have considered here. 

Since 1983 [2] it was believed that the chiral soliton approach allows to describe the 
properties of baryons (the B = 1 sector) with accuracy ~ 30% and better. We have shown, 
that quite good description of the binding energy differences can be obtained also for light 
nuclei B = A < 32, when appropriate rescaling of the constant e is made, providing the 
"nuclear" variant of the model. Taking into account the former results on good qualitative 
description of binding energies of light hypernuclei [13], we conclude that effective field 
theories, including the chiral soliton approach, provide good description of certain basic 
properties of relatively light nuclei. 



We can conclude, that besides the B = 1 version of the model, one can consider the 
"nuclear" version of the model with rescaled value of e which is good for the description of 
nuclear properties for A < 30. Observation of the predicted states can provide additional 
support of the validity of this approach which turned out to be very effective in the under- 
standing of the properties of mesons and baryons. Studies of other properties of nuclei (e.g. 
magnetic moments, radii of different distributions, etc.) is now one of actual problems. 

No stable multineutron nuclides are found within the chiral soliton approach; all 
states are several tens of MeV above the threshold and can decay strongly into final states 
consisting of nucleons. The states with isospin T > B/2 also can be obtained. They can be 
interpreted as states consisting of nucleons and A~ isobars, and some of them can be negative 
- as n...nA~ . Estimates show, however, that their energies are higher than thresholds for 
strong decays. 

One of the evident drawbacks of this approach is that one cannot describe the vari- 
ations of the isotopical content of nuclides in space: in the quantization procedure solitons 
are rotated in the isospace as a whole. In particular, the description of neutron halo on the 
periphery of neutron rich nuclides, established experimentally (see, e.g. [20]), is behind the 
possibilities of our approach in its present form. The studies of nonzero modes quantum 
fluctuations are necessary for this, since the approach we used corresponds to excitation 
only of collective rotational modes. 

Careful study of the results presented in Tables 2—5 and 7—9 allows to reveal the 
following general feature: for relatively light nuclei, A < 20—22, theoretical predictions for 
the difference of binding energies of isotopes are greater than the experimentally observed 
value s of this difference. For A > 22—24 theoretical predictions for this difference are 
smaller than data, as a rule. One could reach better agreement with data by decreasing the 
value of the constant e at smaller A (thus providing the additional increase of dimensions 
of multiskyrmions), and increasing slightly e at larger A. It means, that the relative role 
of nonzero modes in formation of the size of nuclei slightly decreases with increasing A, 
probably, because for large A the size of multiskyrmion is large enough. This observation 
shows that further more detailed studies and refinements would be useful and of interest, 
although they will not change our main conclusions concerning the validity of the chiral 
soliton approach for description of the binding energies of nuclides. 

We are indebted to E. Konobeevsky, A. Reshetin, and Yu. Ryabov for discussions on 
the initial stages of the work, stimulating our interest to the problem considered here. The 
work is supported by RFBR, grant no. 01-02-16615. 
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